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ABSTRACT 

The effect of a leading edge flap on the aerodynamics of a low aspect ratio delta 
wing is studied analytically. The separated flow field about the wing is represented 
by a simple vortex model composed of a conical straight vortex sheet and a con- 
centrated vortex. The analysis is carried out in the crcsa-Sow plane by mapping 
the wing trace, by means of the Schwarc-Chrictoffel transformation into the real 
axis of the transformed plane. Particular attention is given to the influence of an- 
gle of attack and flap deflection angle on lift and drag forces. Both lift and drag- 
decrease with flap deflection, while the lift-to-drag ratio increases. A simple coor- 
dinate transformation is used to obtain a closed form expression for the lift-to-drag 
ratio a a a function of flap deflection. The main effect of leading edge flap deflect ion 
is a partial suppression of the separated flow on the lets ids of the wing. Qualitative 
comparison with experiments is presented, showing agreement in the general trends. 
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perturbation velocity potential 
vortex strength 

exponent in the Schwan-Chrisioffel transformation 
angle of attack 

half apex angle of a main delta wing 
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coordinate of the physical piano 
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coordinate of the transformed plaue 
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tranfomration coefficients 
uniform velocity 
complex potential 

the ratio of the span length of the main v/iag to the total span length 
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normal vector of the 3ap surface 

resultant force 

resultant force on the main wing 

resultant force on the flap surface 

velocity components 

lift coefficient 

drag coefficient 

pressure coefficient 

normal force coefficient 

lift increment ; C L ;i ap - C L without flap 

normal force 
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1. INTRODUCTION 

It is a well known phenomenon that above the lec3ide of a low aspect ratio delta 
wing at incidence there occurs a vortical flow originating from leading edges, and 
that this vortical flow is stable up to certain angle of attack. This vortical flow causes 
an enhancement of lift, called vortex lift, especially needed during maneuvering aa 
well as in take-off and landing. Leading edge flaps have been suggested as a control 
device of this vortical flow to retain the necessary aircraft stability and control. 
This application of the leading edge flap Is suitable for subconic-transonic tactical 
delta wing aircraft, as well as aircraft having both supersonic cruise efficiency and 
high maneuverability. 1 

Coe and Weston 2 made experimental investigations on the effects ofleading edge 
flaps. From their flow visualisation it was found that leading edge deflection was 
effective in partly suppressing the formation of leading edge vortices and promoting 
attached flow condition. Hoffier and Rao 3 conducted experiments and theoretical 
studies using the vortex panel method, including the effect of tabbed flaps. Their 
results indicated an increase of lift-to-drag ratio. BVora the experimental results 
known so far we conclude that during flap deflection the drag 1c reduced considerably 
and the vortical flow is partially suppressed. The former contributes to aa increase 
of lift-to-drag ratio and the latter causes a detrimental effect on lift. The second 
effect has been overcome by the use of flap tabs. 5 ’' 5,5 

Although many theoretical studies on the delta wing configuration have been 
conducted, 6,7,S ' 9 the flapped, delta wing has not been analyzed outside tb.c panel 
method context. 3 Me :cover, the panel method does not provide a satisfactory 
explanation for the experimentally observed changes in the lift-to-drag ratio. In 
the present study we attempt to assess tbs trends of the effect of a leading edge 
flap on the aerodynamic characteristics of a slender delta wing. For mathematical 
convenience the conical geometry of the delta wing and the flap are chosen ns 
shown in Fig. 1. The analysis is conducted in tbs cross-flow plane normal to 
the wing with the assumption that the Sow field is also conical. For this purpose 
the simplest model for the separated flow 9 is chosen. This model is composed of 
a concentrated vortex and a straight vertex sheet, as shown in Fig. 2. In this 
model the analysis is carried out in a plane normal to the wing surface, which is 

l 


* 





shown in Fig. 1. By means of a simple coordinate transformation the aerodynamic 
forces are resolved into lift and drag. A closed-form expression for the lift-to-drag 
ratio is found, which indicates that this ratio must increase with flap deflection. 
Comparison between theory and experiment chows generally cimilar trends for the 
aerodynamic coefficients. 










2. MATHEMATICAL MODEL 


2.1 Methodology 

To study the qualitative trends of the effect of a leadin $ edge flap on the aero- 
dynamic characteristics of a low aspect ratio delta wing, the simplest model for the 
vortex sheet, suggested by Brown and Michael 0 , is chosen. This model, as shown in 
Fig 2, simply assumes the rolled-up ohcar layer to consist of a straight vortex sheet 
and a singular vortex core and requires thin singularity system to be force free in a 
global sence. 

The trace of the wing in the cross-flow plane is transformed into the horizontal 
axis of the transformed plane by the Schv/ars-Christoffel tranformation. The ad- 
vantage of this transformation is that the boundary condition on the wing surface 
is satisfied automatically in the mapping plane. 

In this chapter the following effects of flap deflection are calculated and compared 
with experiments. 3 

© The influence of pressure forcer acting on the flap surface on the iift-to-dvag ratio. 
© The drag on the main wing and the flap surface. 

© The lift increment on the main wing and the Sap surface. 

• The influence of flap deflection angle on the vortex position and the fcrce normal 
to the wing. 


2.2 Analogy between 2-B unsteady flow and 3-D conics, flow 

In this section a comparison between the three dimensional vortex sheet growing 
conically in space and the two dimensional vortex sheet growing in time is made. 
For the detailed mathematical procedure, refer to the paper by Kiichmann and 
Weber.® 

For an incompressible inviscid and iiTotational Sow the governing equation is 

V 2 $ = 0 (1) 

where $ is the perturbation of the velocity potential due to she wing and the 
associated vortex sheets originating at the leading edges. 


i 





lor a conical Cow the perturbation of the velocity potential and the geometry of 
the vortex oheet can be written 

®(W) = z<f>[?,0) (2 ) 

S(x,r,&) = r - z f(0) ( 3 ) 

where r = r/i is the conical parameter. 

In conical coordinates (1} becomes 


IS + (X + + ±2^ „ 

rdr 'fir 2 p2^2 U 


The condition for the vortex sheet to be a streamline is 


V - VS =-. 0 , 


Combining (S) with (2) and (3) gives 


I )' - % + !& - - 


v fir/' fir ‘ / fit' " u 

where - r(d$/dr) y*. 

The Bernoulli equation on the vortex sheet, postulating that no pressure difference 
across the vortex sheet exists, has the form 


For an unsteady two dimensional flow, the potential function and the geometry 
of the vortex eheet can be written 


= y(t) 

S(x, ,3) - r - Jt(±) a m 






.i l.. . : 


where? = (%)/({') is Prandtl’o unsteady similarity variable, and T and R 
are a characteristic time and length. 6 

For the unsteady case we finally have the velocity potential equation as 


W 4. 4. l * 7 * _ n 

?o? + ar 2 r*ao* 


( 10 ) 


Boundary condition for the vortex sheet 


d* f'dfi _ 
m / — — + — — = 0 
J dr f 09 


( 11 ) 


Bernoulli equation for the vortex sheet 


P 

p/2 




Let's compare the velocity potential equation, the boundary condition and the 
Bernoulli equation on the vortex sheet for the usnsteady case ((10),(I1),(12)) with 
those of the conical case ((4),(G),(7}). First, the equations for the velocity po- 
tential in the two cases become identical if r <X X. Next, the boundary con- 
dition and Bernoulli equation on the vortex oheet become equal if tn = 1 and 
tf x = <t> — r(0<t>/dr) <£ 1. Finally we must get the relationship between x and t. This 
relationship can be obtained by returning the potential equation in both cases to 
the original variables. This gives 


/x 2 1 s 1 1 \ 

s ^^2 $ 7 * + ~ < Pr — 2 ~<f>xr *t- <f>rr + ~<Pr + = 0 


(13) 


4 {~~T?4u 

\m-r 


2 — 3m , 12m 2 — 3m +1, t 2 — 2m , _ t 


rm 


Vr + 7t 


~2^06 + + $rr) = 0 . (14) 

r* r • 



With m = 1 and z—Ut these two equations become identical. 

Now we see that the three dimensional conical flow could be calculated from the 
two dimensional unsteady flow if the following criteria were satisfied 

(1) r < 1, meaning that the Cow is geometrically slender. 

(2) Vj <£ 1, meaning that the flow is Gcrcdynzmicalltj slender 

If these two conditions are satisfied, we can celve the conical Sow problem from a 
2-D unsteady problem with the substitution of s = Ut and m = 1. 

The flow around a low aspect ratio delta wing satisfies these two conditions 
except near the vortex core, where v x ia close to U. Then, this analogy can 
be applied to slender delta wings with caution. Moreover, in concentrated vortex 
models, the axial velocity becomes singular at the core location, making the analogy 
more restrictive. 

2.3 Vortex modelling 

In the present model, composed of a straight vortex sheet and a singular vortex 
core as shown in Fig. 2, the vortex sheet and the vortex core do not follow stream- 
lines; that is to say the vortex sheet and the vortex core are not independently force 
free. Thus, Joukowsky forces act both on the vortex sheet and the vortex core. In 
Brown and Michael’s method these two 'forces are required to balance each other 
out, meaning that the whole system is fores free. 9 Uring the unsteady analogy 
derived in 2.2 and from Fig. 2c, we can find a force free condition for the whole 
vortex system. In the unsteady flow, the force free condition is 

p'j - v„„) + - «•„■) = o , (is) 

where F,- !3 the vortex strength, 
cr vi is the position of i l!i vortex, 

V a „ is the velocity at the i th vortex position excluding its own induced velocity, 
a 0 i is the origin of the straight vertex sheet 

The first term is the force acting on the vortex core, the second is that acting 
on the vortex sheet. The term da v ;fdl — is the relative velocity causing the 



Joukowsky force on the concentrated vortex core. The rate of change of the 
vortex strength is d(sI\-)/df and this amount of vorticity i b supplied through the 
straight vortex sheet having uniform strength, which connects the points frc“i a 0 { 
to ff vi . 


EYom 2.2, the relationship between the 3-D conical flow and the 2-D unstead j 
flow is 

A A 

(1C) 


± = u± 
it dz 


and from the geometry of a delta wing, the conical relationship gives 


h = 


se 


(17a) 


— _ 5 = ir 

dx “ x ~ h lj 


(176) 


where £ stands for any quantity associated with the kinematics of the flow field, 
and k is defined in Fig. 3 and e is the half apex angle cf the main delta wing. 


Then the force free condition on the i (!i vortex system can be rewritten 


U e. 

T (2-,i 

i &oi) ~ 

(10o) 


d iT.- 


5 


(18&) 


where F is the total complex potential of the flow field. 

In our case we have >. system of four vorticea fed by vortex sheets emanating 
from the leading edges and from the Sap junctions. Since only the symmetrical 
configurations are treated, the analysis considers one half of the cross flow plane. 
Thus cr^x = * (fc + (1 ~ k)c’ s ) for the leading edge vortex and <r 0 « = * k for the 
flap junction vortex. Then the force free conditions for both vortex systems arc 


))U„ 


d_ 

dcr' 


( 19 ) 
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Uc 


rfavi+ik) 


d_ 

da 


( F -—logic- a v2 ))\ 

\ 4M /I ' 


( 20 ) 


The derivation of the complex potential is facilitated by transforming the trace 
of the wing, together with its line of symmetry, into the real axis of a transformed 
plane. To perform this transformation the Schwarr-Christolfd transformation is 
used. The mapping procedure is shown in Fig. 3. 

The general form of Scwars-Christoffcl transformation is 


da 

dr) 


v - v «) 


a; 


i-l 


( 21 ) 


In integral form 

rt) '* 

a = C + A }T (f? - 7i) fti d V > (23) 

where a,- = kjz — I , /;,• is the angle at the i ih vertex of the polygon, and n; ir. the 
corresponding point of a in the mapping plane. A and C are complex constants.. 
From Fig. 3, we get 

ati = x , &2 = Jr/2 , /; 3 = jt - 5 , hi = 2x , k$ — tt -f 6 , !:q -- rr/2 , (23) 

then 

aji = 0 , 02 — -1/2 , aj = - 6 / rr , 04 = 1 , 05 = V 57 » ~ -1/2 • (2-1) 

From Riemann’s mapping theorem we can arbitrarily choose two points. Let those 
two points be 

174 = 0 (25a) 

Vo - r,2 (2E6) 

3 




h. J 


These conditions imply that the leading edge ia transformed to the origin and the 
initial point of the integration ia taken arbitrarily at tft. Requiring that the velocity 
in both planes be equal at infinity, we have 


dff 

df} ff-*oo 



(26) 


With (25) and (26), the transformation function becomes 


<j 


P 7 

ht 72 )(7 - 7s) 



(27) 


Theconst?-ntc i 72 » 7s> *7s» 76 can be obtained by matching the four corresponding 
points in the physical and the transformed plana 


cr s = tfc 


73 

(28s) 

<r 4 « ik{l + (l- k)e iS ) 


0 

( 28 h) 

os = ik 

— ♦ 

75 

(28c) 

0 

ti 

b 

— *> 

76 • 

(23d) 
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Here lengths are normalised with the half spun length a(x). 
Applying (28), 


k+ 

t 

P* « 

/75-f 

\6/* 

) « 

In* vW - OW ~ 72 ) 

'73 — 

k+ 

e ( 

f >75~e\ 

<V* 

cu 

f! 

h* n/(7o-OU*-72) ' 

^“73' 

k- 

P » 

nis- C 

,6/r 
) *<■ 

« 

/o \/(70 ”£)(£ - 72) 

H *“ 73' 

it- 

€ 

/£ — 75' 

,5/r 
) df 

« 

4s \/(76 “ £)(£ ~ 7z) 

W-73. 
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(29a) 

(20&) 


(29c) 


(29,1) 





Thi 3 system- is solved numerically for 32, 7 s> *?s m ohown in Appendix 1 . 

FVom Fig. 3 we nee that the complex potential in the transformed plane is the 
sum of the potentials of a unform velocity and that of four vortices, two of which 
are the image vortex systems needed to satisfy the boundary condition on the wing 
surface. Hera we neglect the influence of the vortex 3 heet on the flow field. 

Then the complex potential in mapping plane is 


F(jj) = Ucxt) 4- 



( 30 ) 


and the complex velocity in the physical plane is 


dF_dFd 1 
dcr drj dir 

r • . , * ( 3i ) 

= lua + —■{ — ) - ^ ( — — - — }] 

L 2«r \jj — ij v i tj-r) vi S 2k- Vf? - v } v 2 t? - j &r 


where 

df? _ v/(Y~- f?a)(7~;2 ) ;' »7 - Vi \ s f* f32 » 

do- 9 '17-95^ * ' 

According to ( 32 ) the velocity becomes singular at the leading edge, ij = 0 , and 
at the flap junction, rj = 95. To avoid this singularity wc apply the Kutta condition 
at these two points. The two equations for the Kutta condition are 


Ua (— — ) - ^ (~ L) - 0 , 

2ir Vr7t l 9rl' 2?r \rf v2 W 

Ua + f - ?^r-) - — ( - ~) = 0 

2?c \i 75 -n..i 2x \n 5 - n.. 2 rj 5 - 


(33) 

(31) 


With the complex potential F the force free conditions given by (IP) and ( 20 ) 
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become 


U £ 
k 


(2yrf+*(* + (i - /; ) c " 15 )) = 


— j 
k ' 



r , »Ti 
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»t 2 / 

1 

— 1 )1 


c/cs — — 

2jt *7t,i 

-»7vi 

2 ir V 

*7vl ~ 

*7f2 f lvl ~ '7t’2/ J 


Z' \A»7ei ~ 

• ^cK^l 
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- »/fl) “ 

log(v * 

- Cfl) 
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.U ,T *. 1 1 


[ 0 2w \tiv2 - 

«7v2 

- s 7t>x 

/ 2sr t/y'j - ‘7 u2. 


A/fai*’ 


2 ~ fJ2) 

! ntvv. 

- %\ 6 /*\ 


l 
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— *75^ / 
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log(cr 

~ <?&) 

] 

J *y— 
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( 35 ) 


Combining (33) and (34) with (35) and (36), vre get the following two complex 
equations for the vortex positions. 


e +<(*+(»-*)«"“)) = 

r, _ P-g ■ _ C -' 1 ( i X (37) 

[ CB - AD TM - rfv 1 ~ ^ '‘‘M “ Tt-2 *?vl “ 






c 


ka 


(2cr v 2 + •*) = 


1 + 


D-B 


(_■ 4-) 

\n,. < > - n.,i fJuS — Aui' 


+ 


C-A 


1 


©L 


CB- AD \ rj„. 2 - <7t-i f/t-2 - Vvl ’ CB- AD t) v2 - T} v2 

A -C d _ nv7 ) _ / og (or - c t2 )) 


(38) 




CB~ AD dff' 




where 

A- ^ (3Sc) 

Avl Ai>l 

J3 = J L (396) 

J 1v2 7 !v2 

c = - L __ (33c) 

% - *7wl A5 " At* 

£ = A . (3Gd) 

‘75 “ Ar2 AS ~ Ac2 

In the absence of flaps the only governing parameter would be the ratio of e/ft . 
In the present case two additional parameters ,5 and k, are included in the group 
£/(/.* a). The flap deflection angle 8 does not appear explicitly but through the 
transformation parameters 1721 A3 > AS a ^d A&> which are functions of 0 and k. 

A Newton-Raphson meth d is used to solve the ; beve two equations, the detailed 
numerical procedure is given in Appendix 2. 

With the vortex positions known we get the vortex strengths from 



ill ff-rJL ua 

2jt CB - AD 


pf. 


' b 

r i 


i r 2 C — A „ 

17 " CB - AD ' 


2.4 Ac-rodynamic forces 



(40a) 


(406) 
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2.4.1 Rcenltfint foyco 

To calculate the resultant force on each surface we consider the pressure coeffi- 
cient. 

C p = l + Q 2 -(u/£f) 2 -(«/17) 2 -(«;/^) 2 , (41) 

where the v and u» components are given by 



= Re (i V) ( 43c ) 

here a = y + i z and V = ( dF/dr)){dr)/da ) i.e. the complex velocity in the 
crosa*3ow plane. 

The u component is beet obtained from the conical flow condition 
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The presence of the feeding sheet induces a discontinuity in the u component of 
the velocity, aa demonstrated in Appendix 3. 

The resultant force is given by 

R = f p da (45) 

Jeach stir fact 


Thi3 integration is carried out numerically in the physical piaas. 

2.4.2 Lift-to*Drag ratio 

Most experiments 2 ^’ 5 show an increase of iift-to-drag ratio with Gap deflection. 
Before we calculate the lift-to-drag ratio we must consider the normal vector on the 
flap surface. The rotation of the leading edge flap about the flap hinge axis causes 
the normal vector to have a component opposite to the free stream. Associated with 
this component there is a force tending to decrease the drag, which is responsible 
for the increase in CifC& . This component can be obtained by performing tie 
following transformations: 

Let’s take the longitudinal axis of the wing aa s, the axis normal to the wing 
surface as y and the starboard direction as ~z, as shown in Fig. 1. e x ,e y ,e a 
denote the unit vectors along axis s t y, z . With the rotations shown in Fig. 4, 
these vectors are expressed as 


/ex, \ / cost 0 jin A / e x \ 

K. « ? i 0 M 

\ c a , j \-jih£ 0 cost) \e 3 / 


Identifying n with the normal vector to the flap surface v/e have 


n = e 


ya 


(46) 


/°x* \ / 1 0 0 \ / °x i \ 

I e ya | = 0 cos S sm <5 J I *i y , ) . (4?) 

\ e c* ) VO —sin 5 cos 6 } \e* t / 


(43) 
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Then 

a = cos 6 Cy l + sin 6 e B , 

= cos S Qy t + sin 6 (sin c e* + cos e e a ) (49) 

~ -e sin 5 -!* cos 5 o y -i- oin 8 c B 

The first term in the last line shows the forward component of the normal vector. 
This term contributes to the increase of the lift-to-drag ratio when the flap is de- 
flected. The first component, —esinS, gives the thrust force and is of higher order 
than the other two. As the flap is deflected the vorta: coscto closer to the sving 
surface, causing a considerable change in drag. 

On the main wing the resultant force is normal to the wing surface. Then the 
lift and the drag components are 

L w = R w (50s) 

J?iy = i?cy (i (c*C5) 

where R w is the resultant fores on the main iving. 

On the flap surface there are all three components of force, as indicated by («D) 


■? 

i! 

1 

n\ 

Ci 

c* 

(51a) 

Rpy — Rp cosS 

(516) 

Rp z = Rp e oinS 

(Sic) 

where Rp is the resultant force on the Sap surface. 

Since the Rp. component ia cancelled by symmetry, only Rp 2 
to lift and drag 

and Rp v contribute 

Lp — Rpy — Rp x a 

(52a) 

Dp = Rpy a + Rp x co-iS 

(526) 


FYom (51a) and (51b) 

Rp x — — « tan5 Rp v . (53) 
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From (51) ,(52) and (53) we get; the drag and lift on the flap 

Lp = Rpy (l + «a ianS) 


Dp = a Rpy f 1 — -- tanS^ (546) 

The total lift and drag are 

L — R w + Rfu(1 + £ ° foni) (55) 

D — a + Rpy(j- — — ians)^ . (56) 

(55) and (56) show the influence of the forward component of force on the lift and 
the drag. Due to this thrust component the lift increases by an amount of order 
cr 2 , whiie the drag ia reduced by almost the came order as the drag itself. This 
increases the lift-tc-drag ratio during flap deflection. 

For moderate deflection angle the lift-to-drag ratio is 


-~tan S ^ + 0(e 2 ) . 
.or / 


a|y Jr Rv + Rpy 


Numerical results show that the lift itself also decreases with flap deflection while 
the ,lift-to-d~,g ratio increases. 

2.4.3 Normal force 

As mentioned in the previous section the lift decreases with flap deflection. By 
considering the normal force we can infer the reason for this reduction. 

The normal force can be calculated by considering the momentum change in the 
cross-flow plane 




Applying Stokes’ theorem 


N- -pU j - 


1C 



(60) 


In terms of the complex potential and the real stream function 

N = - pU f fon (F da ~i § da) 

Je 

Since the contour is a streamline 

//= -pU jflm (F da) 


( 61 ) 


Substituting for F 


j»— + fc,(^)) £* (62) 


To perform this integral we first express the part containing the logarithmic terms 
on a circle of infinite radius, where da/dq = 1 


J c \2tt \q-q v i / 2jt \q ~q v ^} dq 

= f (p.iUl=*±)-.p h J!L rJM))*, . 

/r=co\2*r \q-t lvl J 2s: J \q~q v 2 ^j 


( 63 ) 


We see that in the transformed plane the normal force arising from the logarithmic 
terms does not depend explicitly on the wing trace, but implicitly through the 
vortices positions and strengths. This equation shewn that the normal force is just 
the sum of the rate change of impulse of two pairs of vortices. Then the result could 
be simply obtained by considering the impulse of a pair vortices, or by integrating 
around the vortex sheet aa shown in Fig. Ga. Following this icond approach give-3 

/ ~ l ° 3 ^ ” lfv ^) dn 
sTi r f n ' 1 „ 

= 2^r [y_. l °3i r > ~ »7«l)lp”^(i7 - qvl )!p-» <*n 


( 64 ) 


rn,i _ 

+ / log(q - «7ti)lp ~ log(q - *? vl )| p dq 

J >ht 

0?t>l " *7;>l) ■ 
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-i -t : i . 

Then the normal force originating from the logarithmic terms is 

N c - Im pU{ Ti(rj v i — tj 0 {) ~ r 2 (s?f-2 “ »7t;2)] • 




(65) 


This force i3 associated with the separated Cow field about the wing. 

The component of normal force from the the first term in (G2) has a singularity 
at infinity. This component could be identified with the force produced by the wing 
if the flovr were attached. Integrating dong the wing trace, as shown in Fig. 5b, 
we find 


r *7 — *75 \ S I* 


f 5 ( '1 ,<0 ' * j- 

Jc VUT- *70 )(*7 ~ V2r^ 

= P I /r?-f75\^ g | d 
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where subscript p represents the principal value that means all values must be 
calcualted in the same Riemann surface. 



Then the total normal force is 


// = im (pu r £ (i 7 ,. x - sj vl ) - pU r 2(^2 - 

— pU* a Isa [ — ^- 2 f 2-^5)% 

•/c %/(9 " »?$)(*? -.>72) '*? ” f /S' 


(67) 


The normal force is composed of two parts; the normal force from equivalent at- 
tached now and that from the separated Cow. 

Dividing this total normal force by the dynamic premure pU^/% and the pro- 
jected area k a(s) (A + (1 — k)cos 5)1 z , we get the normal force coefficient 


= r 


b k + (1 — k)eos5 

I® 2Fi(n,,i - ?/ t . 1 ) — 2F 2 (q, 3 — Tf v 2 ) — 


* (cc + /: + d ,r 


>< 7s ' ~ — * 72 ) ‘*7 


(“>" /c 41 

'•> D / J 


(60) 

where tic- vortex intensity and q have been normalised with U a(s) and a(s) 
respectively. 


Here ail the values of the integrand are the principal values. This shows explicitly 
the effect of the Cap angle on the normal force. The region from 173 to V 75 is the 
flap surface. 

For flap deflection angles much larger than the angle cf attach, it is conceivable 
that the vortex may be located underneath the wing surface . 6 In this case the 
dominant component in the cross flow plane is not U a but the projection of U 
on a plane normal to the flap surface. This fact imposes a limit to the flap angles 
for which this theory is expected to apply. 
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3. SS3ULT3 jawa COKCLUSEOWS 

3.1 Remits 

3.1.1 CoaforsauS aaeppiBg 

The following tabic chov?a tie SracnsKTEStion cotScicaia -»j2» »i3« *J4 r>ndt 2 $ for 
■varying Hap deflection angle 5 , sitin';-— 0.5 aad e •= 35 - 95 ° . 


£(deg.) 

% 

^3 

^75 

% 

40. 

-0.6369 

—0.5268 

0.9785 

1.301 

36. 

— <5.6723 

-0.5441 

0.0650 

1.277 

32. 

-0.7051 

-0.5722 

0.9533 

1.251 

28. 

-0.7455 

—0.6004 

0.9355 

¥ 

24. 

-0.7024 

-0.6204 

as* S3 

2.105 

20. 

— 0JUC1 

— CuSJraS 

00357 

13SS 

16. 

—D.B5S3 

—©.6533 

0JB7S6 

2.234 

12. 

. — ojbbsa 

—0.7205 


■H « ir.-r*. 

L.&UA 

S. 

— OJ3'i3S 

-0.7371 


2.053 

4. 

—£115544 

-D.7029 

OS22& 

2.035 

2. 

-Q.9322 

-0.7755 

06020 

3.028 



Table I 




We sea that as 0 — » G . V72 — * — 1 aad * 1 , v?.bkh ccarsaj and io She Brovra 
and Michael case. 

3.1.2 AweodyTuaadE fcrcca 

To compare with the experimental ie3ultu of ref. }3] a wiug with c — 15.33° and 
& -•= £}.G is chosen. The angle of attach h varied from 0° to 20°. 

Tig. 6 eiicwo the contribution to the drag of the main win^. B css be oean that 
the model caeraitimatefl the drag, particularly at higher angles of attach. 

Fig. 7 chows the contribution to the drag of the flap. Notice that this component 
of drag can roach negative values. Thin ia due to the forward force component 
discussed above. 
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Fig. 8 and 9 show the increment of lift on the main wing and flap respectively. 
The abrupt departure from experiments observed at about 12° appears to suggest 
vortex breakdown in the experiments. 

Fig. 10 shows the evolution of the resultant drag. The negative values of drag 
arise from the conical assumptions, where no trailing vortices and associated drag 
are considered. 

Fig. 11 depicts the total lift change. Again, the sudden departure from the 
experimental trends may indicate that those measurements were affected by vortex 
breakdown. 


Fig. 12 shows the change in the CijC p ratio. No experimental values were 
available. The no-deflection case b included for comparison. 

Fig. 13 shows the positions of vortices for varying flap deflection. For flap 
deflections less than 18° the secondary vertex almost disappears. The vortex 
positions in the physical plane are ehowu in Fig. 14. We notice that flap deflection 
causes the secondary vortex to move away from the surface and the leading edge 
vortex to approach the surface. 

Fig. IS describes the strength of both vortices aa function of flap deflection. Flap 
deflection strengthens the secondary vortex and weakens* tho. leading edge vortex. 
The strength of leading edge vortex reduces by a half in 40 3 flap deflection. 


Fig. 16 shows the different components of normal force. We notice that flap 
deflection acts primarily on the leading edge vortex component of the normal force. 
The secondary vortex has a negligible effect the total force. 


3.2 Conclusions 

A conical delta v/iug with leading edge flap3 wap analysed using a simple vertex- 
feeding-sheet model consisting of two pain of vortice-3 with corresponding vortex 
sheets. One pair originates at the leading edges and the other at the flap junctions. 
The leading edge vortex pair is much stronger and responsible for most of the lift and 
drag changes that take place during flap deflection. In agreement with experimental 
evidence, it was found that flap deflection causes a decrease of lift and drag, and an 
increase of the lift-to-drag ratio. A closed form expression exhibiting the lift-to-drag 
ratio increase was derived. The rapid reduction in drag is attributable to a thrust 


force component developed by the flaps, while the main reason for the lift reduction 
is the partial suppression of the vortical Sow. There ia agreement with experiments 
in the general trends. 

The two main sources of limitation of the theory are the c redness of the vortex 
system representation, and the assumption that a cross- Dow plane can be suitably 
defined. The latter aspect becomes relevant in highly warped wings, which in this 
case amounts to flap deflections that are large and comparable to the angle of attack, 
where the notion of a plane normal to the wing becomes meaningless. 
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APPENDIX 

A 1. Numerical procedure fo? conformal mapping 

A 1.1 Welfo method 

',Volfe : o method is just a variation on the secant method for systems of nonlinear 
equations. This method has the advantage that as the solution becomes closer to the 
root, the convergence speed is much faster than that of NewtamRaphcoa methods. 
Fur a single equation the secant method is equvalent to solving the following two 
equations, with two initial gucssea 


Pl + P2 = l 

(41.1a) 

Pl/(*n) + F2 “ 0 

(,41.16) 

The solution is updated as 


S n+X ~ Pl s f»~i *h P2 T -n 

(Al.lc) 

For a system of n equations a «e» of ft 4- 1 equations must be solved with n 4 1 

vectors of initial guesses, 


JM? 

ii 

( 41.2a ) 

n+i 

Spj fi{si J >X2 J t'",x n >) = 0 t = 1,2, -.a . 

;= i 

(41.26) 


For the solution update we ignore the vector x J which gives the maximum residue. 


n+1 

*tep = P/ *frp , (41.3) 

j-l 


where x* = (ri J , x ? , • • • , x n J ). 
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This procedure is repeated until the maximum residua becomes smaller than a given 
value. 

A 1.2 Integration 6ft oinguior points 

The integrand of the mapping function is singular at the points £ — tj 2 , tj St q s . 
To perform the numerical integration, the integrand is expanded to 1st order in 
Taylor series about those points. For example, at £ ~ »J2 


f'** 1 e ( fis-t yi* dfi 

= r 

Jo 


f + *?2 


/ *?s -- f — g;a \ c / y . 


S- *12 / t?5 ~ *12 \ 6 I* /* J_ 

>/«76 ~ »72 ^3 “ m' Jo VS 

= ggs nm-th y/ g ^ 

Vh<s - m W3 -%•' v 


d f 


where f = { - vj 2 - 

Romberg integration is used on the rest of the range. 


A 2. N‘fui.:e?ka! procedure for ika vertes: positions 
A 2.1 Newton - Rnphscm saeibod for complex fxmefciosss 
Consider the system of real equations 


£(::) = 0 

The general form of the Newton-Rapheon method is 

Xn+l = ZZ n - £(»».) 

where J[f(x n )] is the Jacobian of £(?:) at x n . 

In the case of a system of complex equations 

X = (xi +m, *2 + * 1*2, •••,=( 1 -MVn) 
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(A 1.4) 
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(A. 2.2) 


(A. 2.3) 
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Applying the scheme to the real and imaginary parts of ? we get 
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(A.2.1) 


(A.2.5) 


matrix. The reason why the Nev/ ton - R&phsoa method is used here is that for Wolfe 
methods we need one more initial guess than for the Newton - Raphaoa method, 
which givta 6x6 matrix. A more important reason is that the Wolfe method needs 
very accurate initial guesses for good convergence. 


A 2.2 Integration path for the mapping function 

In using Newton - Raphcon method to solve the equations for the vortex positions 
we need to calculate the corresponding points between the physical plane and the 
transformed plane by carrying out the transformation 


c 


[ n •n 
hi VTn - j? 2)(»? - tie) 



(A2.6) 


If we take the integration oath as a straight line connecting J 7 2 with the assumed 
vortex position, we encounter the following problem: when the vortex position is 
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very close to the real axis, which b the case when a b small, the integration 
path almost passes through the singular points ;,>$ and ij<j. To avoid this possible 
numerical singularity, the integration path is taken in the following way: 

t)z — * Tjz + » 2 — » assumed Re f] v + * 2 -» assumed v;„ (.,12.7) 

The maximum imaginary value is taken as 2, which ir; large enough to calculate the 
vortex positions. 

A 2.3 Velocity cfc fclia vortex location 

(37) and (30) in chapter 2 contain terms of the form 


dj 

dcr 


[M? - rj vi ) - log(<T - or v ;)]^_ 


(A2Ma) 


which can be v/ritten as 

_d_f 

dq 


~ l°9(° " ^vi)l j~| 


(A2M) 


In terms of Snite differences this quantity can be approximated by 



r?,.,- 4- Ar?,,,- - r lv{ 
a{n v ; + &«?„•) - cr vi 


[on tyfi "7 vi 

9 cr[>ivi ~ Aq.,.) - a ui 



[\ cr(j? t .,- + Ajjn) -o vi n \ 


ds >u< 


(>12.9) 


This is done in the transformed plane to avoid inverting (A 2.6). 

A 3. Derivation oi* the axlc.1 velocity component 

In 2.4 w« derived the axial velocity component under conical (low assumptions 



(A3.1a) 




• yv 


au»](A3.1t) 



Here we have logarithmic branch points at i]vl> ?jv2> t/ vll lj v 2 - To numerically 
evaluate functions across branch cuts care must be exercised since computers usually 
give the argument of a complex variable between — * and if while the value of interest 
here is between 0 and 2s\ Ae an example we consider two branch points at and 
>f vl , aa shown in Fig. 5a 


login -nvi) = | 


login - nvi)\p— login - t7ci)i c +i 2 s- 
login - *7fi)lp+» 2s- = login - s?i,i)j c +s ‘In- 


in i,H,m,rv 
in V 


(A3.2a) 


login - Vvi) = | 


login - 57 vi)! P = login - nvdU+i 2{r 
login - «7vi)lp+» 2a- = login - 7vi)l c +i 


in I 

in n.m.iv.v 


(A3.26) 

where the subscript p means the principal value and c means the value calculated 
from the computer. Actually the values are calculated only in region IV, V and half 
of III due to the symmetry of the flow. 
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A 4. Fortran program fop tho ScLayari^Clii-'iataSai troj cf&rmatlon 



C700 

1 



IMPLICIT REAL*8(A-H.0-Z) 

DIMENSION SS(10, 1C) .WORK (10) ,0A(10. 10) . Z (10) , IPVT(IO) ,XF(iO) 
,FNRM(10) ( X(10.10) ,XZ(10) .5(10) ,XX(10) , VX (10) , SX(10) 
DATA IMAX , JMAX/4 . 5/ 

DATA PI/3.141592654/ 

DATA ERR, EPS/1 .D~8, 1 .D-2/ 

DATA XRTO / .6153846154 / 

WRITE (6 , 700) XRTO 

FORMAT (5X, 1 FLAP ANCLE = \D15.6) 

DO 1 J=1 , JMAX 

READ (35. *) (X (I , J) . 1=1. IMAX) 

WRITE (30, 121) 

READ (5, *) DLTAD 
DLTA=DLTAD * P 1/180 . 

N=0 

DO 2 J=1 , JMAX 
DO 3 1=1. IMAX 
VX(I)=X(I.J) 

CALL MATEL ( VX . IMAX . JMAX .PI, XRTO . DLTA . ERR . EPS . SX) 

DO 4 1=1. IMAX 
SS(I.J)=SX(I) 

SS (JMAX. J) =1 . 

CONTINUE 
DO 6 1=1. IMAX 
Q (I) =0. 

B (JMAX) =1 

CALL DECCMP ( 10 . JMAX . SS , COND . IPVT. WORK. OA, Z . ZNQRM) 

CALL SOLVE (10. JHAX.SS.B. IPVT, XZ) 

DO 7 1=1, IMAX 
XX(I)=0. 

DO 7 J=l, JMAX 
XX(I)=XX(I)*X(I,J) *XZ(J) 

CALL NORM (OA. IMAX, JMAX. FMX.FNRM. XLOG. JC) 

IF (XLOG . LT . - 10 . ) CO TO 0 
IF (N.GT.SO) CO TO 9 
DO 10 1=1, IMAX 
X(I.JC) =XX(I) 

N=N+1 

WRITE (6, 500) N, (X(I.JC) .1=1. IMAX) 

X(l, JC)=-DABS(X(1,JC)) 

X (2 , JC) =-DABS (X (2 , JC) ) 

X (3, JC) =Di\ES (X (3. JC) ) 

X (4, JC) =DABS (X (4, JC) ) 

IF(X(2,JC) .LT.X(3,JC)) CO TO 300 
CXX=X(2. JC) 

X (2, JC) =X (3, JC) 


i 
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t T. ;'{ ’■ ' 

o . 

ft ..d 




X (3, JC) =CXX 
GO TO 300 

300 IF(X(3, JC) .LT.X(4, JC) ) GO TO 400 

CYY=X(4.ja 
X(4,JC)=X(3.JC) 

X(3.JC)=CYY 
GO TO 400 
400 GO TO 5 

9 WRITE (6. 100) 

GO TO 11 

8 WRITE (6, 500) N, (XX (I) ,1=1. IMAX) 

GO TO 11 

11 WRITE (30. 202) DLTAD, (XX (I) ,1=1. IMAX) 

IF (XRTO.CT. 2 . ) GO TO 130 
GO TO 111 
130 STOP 

121 FORMAT (3X, 'DLTA',9X, 'A2’,14X, 'A3'.14X, 'A5\14X, ’A6',1X/) 

100 FORMAT (IX. ’MORE THAN SO ITERATION') 

202 FORMAT (IX , F8 . 5 , 4 (IX . D15 . 6) ) 

500 FORMAT (IX, 12,4 (1X.D15 .6) ) 

END 
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ft 


ft 
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ft 


ft 
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SUBROUTINE MATEL (X, IMAX, JMAX, PI, XRTO, DLTA, ERR. EPS. Y) 
IMPLICIT REALMS (A-H.0-2) 

DIMENSION X(10) ,Y(10) 

A2=X (1) 

A3=X (2) 

A5=X ( 3) 

A6=X (4) 

A2P=A2->EPS 
A3M=A3-EPS 
A2P=A3+£PS 
A6M=A6 - EPS 


CALL ROMBRG ( 1 , A2 , A3 , AS , A6 , A2P , A3M. ERR . DLTA, PI . XINT1) 

CALL R0MERG(2.A2.A3.A5,A6,A3P,0. .ERR. DLTA. PJ.XINT2) 

CALL ROMBRG ( 2 , A2 , A3 , A5 , A6 , 0 . . A5 , ERR , DLTA , P I , XINT3) ! 

CALL ROMBRG (3, A2 , A3, A5, A6 , A5, A6M, ERR, DLTA. PI . XINT4) I 

Y (1) =XRT0-»XINTl + 2 . fi A2* ( (A5-A2) / (A3-A2) ) * (DLTA/PI) ! 

*DSQRT(EPS/ (A6-A2) ) +A3* (AS-A3) ** (DLTA/FI) ! 

•EPS* * (1 . -DLTA/PI) / ( (1 ■ -DLTA/PI) *DSQRT ( | 

(A6-A3) * (A3-A2) ) ) i 

Y (2) =1 . -XRTO+ 

XINT2+A3* (AS- A3) ** (DLTA/PI) *EPS** (1 . -DLTA/PI)/ 

( (1 . -DLTA/PI) *DSORT( (A6-A3) * (A3-A2) ) ) 

Y (3) =1 . -XRT0-XINT3 j 

Y(4)=XRT0-XINT4- j 

2 . *A6* ( (A6-A5) / (A6-A3) ) * * (DLTA/PI) *DSORT (EPS/ (A6-A2* ) 
RETURN j 

END • j 





SUBROUTINE NORM (X . IM , JM .XMAX , X1SM . XLCG , N) 

ILLICIT REALMS (A-H,0-Z) 

DIMENSION X(10,10) .XKRM(IO) 

XMAX=0. 

N=1 

DO 2 J=l. JM 
XNRM(J) =0. 

DO 1 1=1. IM 

1 XNEM (J) =XNRM (J) >X(I, J) * 6 2 
Ic (XMAX.GT.XNRM(J)) CO TO 2 
XMAX=XNRM ( J) 

N=J 

2 CONTINUE 
XLGG=DLOC (XMAX) 

RETURN 

END 

SUBROUTINE DECOMP (NDIM.N, A, CON’D, IPVT, WORK. CA, 2 , ZNOR.M) 
IMPLICIT REALMS (A-H.O-Z) 

DIMENSION A (NDIM.N) .WORK(N) .OA (NDIM.N) , Z (N) , IPVT (N) 
INTEGER NM1.I.J.K.KP1.KB.XM1.M 
DO 1 1=1. N 
DO 1 J=1.N 
UA(I . J) =A(I , J) 

1 CONTINUE 

IPVT (N) -1 

XF(N.EQ.l) GO TO 80 

NM1=N-1 

AN0RM=O. 

DO 10 J=1.N 
T=0. 

DO 5 1=1, N 
T=T+ABS (A ( I , J) ) 

5 CONTINUE 

IF(T.CT.ANORM) ANORM=T 
10 CONTINUE 

DO 35 K=1.NM1 

KP1=K+1 

M=K 

DO 15 I=KP1,N 

IF (AES (A (I , K) ) GT. A3S (A(M.K) ) ) H=I 
15 CONTINUE 

IPVT (K) =M 

IF (M • NE . K) IPVT (N) =- IPVT (N) 

T=A(M,K) 

A(M,K)=A(K.K) 

A(K. K) =T 

IF (T.EQ.O. ) GO TO 35 
DO 20 I=KP1.N 
A ( I , K) =- A (I , K) /T 
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CONTINUE 
DO 30 JaKPl.N 
T=A(H.J) 

A(M, J) =A(K, J) 

A(K.J)=T 

IF(T. EQ.O.) GO TO 30 
DO 25 IsKPl.N 
A(I.J)=A(I.J)*A(I.K)*T 
CONTINUE 


30 CONTINUE 

35 CONTINUE 

DO 50 K*1.N 
T=0. 

IF(K.EQ.l) GO TO <15 
KM1=K-1 
DO 40 1=1. KM1 
T=T+A ( I . K) * WORK (I) 

40 CONTINUE 

45 EK=1 . 

IF (T.LT.O.) EK=-1. 

IF (A(K, K) . EQ . Q . ) GO TO 90 
WORK (K) =- (EK+T) /A (K , X) 

50 CONTINUE 

DO 60 KB=1 , NM1 

K--M-KB 

T=0. 

KP1=K+1 
DO 55 I=KP1.H 
T=T+A(I,K) -WORK (K) 

55 CONTINUE 

WORK (K) -T 
M=IPVT (K) 

IF (M EQ . K) GO TO 60 
T=WORK (M) 

WORK (M) =KORK (K) 

WORK (K) =T 
60 CONTII'IUE 

YN0RM=O . 

DO 65 1=1. N 

WORM=YNORM+ABS (WORK (I) ) 

65 CONTINUE 

CALI. SOLVE (NDIM.N, A. WORK, IPVT. Z) 
ZNCEM=0 . 

DO 70 1=1. N 
ZNORM=ZNORM-ABS (2 (I) ) 

70 CONTINUE 

COND=ANORM * ZNGRM/’: tfNORM 
■ IF (COND.LT'.l.) C0ND=1. 

RETURN 
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C0ND=1 . 

IF (A(l, 1) .NE .0. ) RETURN 

C0ND=1.E*32 

RETURN 

END 
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SUBROUTINE SOLVE (NDIM.N.A,B, IPVT.X) 
IMPLICIT REAL*6(A-H,0-Z) 

DIMENSION IPVT (N) . A (NDIM, N) , B (N) . X (N) 

DO 1 1=1. N 
X ( I ) =B ( I ) 

1 CONTINUE 

IF(N.EQ.l) CO TO 50 

NM1=N-1 

DO 20 K=1 . NM1 

KP1=IC+1 

M=IPVT(K) 

T=X (M) 

X (M) =X (K) 

X (K) =T 

DO 10 I--KP1.N 
X (I) =X (I) + A(I , K) ft 7 
10 CONTINUE 

20 CONTINUE 

DO 40 K3=l , NM1 

KM1=N-K3 

K=KM1+1 

X(K)=X(K)/A(K,K) 

T=-X (K) 

DO 30 1=1, KM1 
X (I) =X (I) +A(I , K) ft T 
30 CONTINUE 

40 CONTINUE 

50 X (1) =X (1) /A(l , 1) 

RETURN 

END 


SUBROUTINE ROMBRG (N , A2 . A3 , AS , A6 . A , B . ERR . DLTA . PI , RES) 
IMPLICIT REAL'S (A-H.O-Z) 

DIMENSION 2(10.10) 

1=1 

DEL=B-A 

IF (N-2) 30,40,50 

30 2(1, 1) = .5*DEL* (FUN1 (A, A2 . A3 , A5 , A6 . DLTA. PI) 

* ' +FUN1 (B , A?. , A3 . AS , A6 . DLTA , ? I ) ) 

GO TO 10 





Z(1,1)=.5«DEL* (?UN2 (A, A2 , A3 , A5 , A6 , DLTA, PI ) 

+FUN2 (3, A2. A3.AS.A6.DLTA.PI) ) 

GO TO 10 

2(1, 1) = .5*DEL* (FUN 3 (A, A2 , A3 , AS , A6 . DLTA, PI } 

+FUN3 (B, A2. A3, A5,A6,DLTA.PI) ) 

GO TO 10 

J=IFIX(2 . ** (I - 1) ) 

DEL=DEL/2 . 

1 = 1+1 

Z(I,1) = .S*Z(I-1.1) 

DO 1 IC=1,J 

X=A+ (2 . * FLOAT (K) -1 .) “DEL 
IF (N-2> 60,70,80 

Z (I, 1) =Z (I , 1) +DEL ft FUNl (X.A2, A3, A5, A6.DLTA.PI) 

GO TO 1 

Z(I,1)=Z(I.1)+DEL*FUN2 (X,A2.A3, AS. A6. DLTA, PI) 

CO TO 1 

Z(1. 1)=Z(1 , 1) +DEL*FUN3 (X, A2.A3, A5,A6,DLTA.PI) 

GO 'TO 1 
CONTINUE 
DO 2 K=2,I 

Z(I,K} = (4.**(K-1)*Z(I,K-1)-Z(I-1.IC-1J)/(4.**(K-1)-1.) 
CONTINUE 

DIFF-'DABS (Z (1. 1} -2 (1,1-1)} 

IF(DITF.LT.ERR) GO TO 20 
GO TO 10 
RES=Z (1,1) 

RETURN 

END 


FUNCTION FUN1 (X. A2 , A3 . A5, A6 . DLTA, PI) 

IMPLICIT REAL»8(A-H,0-Z) 

FUN1=X* I (A5-X) / ( .3-X) ) 44 (DLTA/P I ) /DSQRT ( (A6 -X) 4 (X-A2) ) 

RETURN 

END 


FUNCTION FUN 2 (X, A2.A3, AS. AS, DLTA. PI) 

IMPLICIT REAL*3 (A-H, 0-Z) 

FUN2=X* ( (A3-X) / (X-A3) ) * * (DLTA/P I ) /DSQRT ( (A6-X) 4 (X*A2) ) 

RETURN 

END 


FUNCTION FUN 3 (X . A2 , A3 , A5 , A6 , DLTA, PI ) 

IMPLICIT REAL 4 8 (A-H, 0-Z) 

FUN3=X 4 ( (X-A5) / ^X-A3) ) * 4 (DLTA/P I ) /DSQRT ( (A6-X) * (X-A2) ) 

RETURN 

END 
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( C 

I c DELTA WING WITH FLAP PROBLEM 

| C WITH ONLY ONE VORTEX FROM LEADING EDGE 

f C USING NEVfTOM-RAPHSON METHODS 

! C WITH DIFFERENT INTEGRATING PATH 

I C 

Q*****#ftAG*i**fcft ft ft*ft ft ***»*>£*<; A* 




IMPLICIT INTEGER (I -N) 

COMPLEX* 16 XI.X.Z.VI 

REAL fc 6 ERR . P I . ALP , AFX . XXI . XK2 , DLT . A2 . A3 . A5 . A6 , CAM , ALPD 

* CL3, CLW , CDW , CLF , CDF , CLT. CDT, CLL , TCL , TCD.RTOR , RTO 

DIMENSION Z(10) ,W(10) .GAM(IO) 

COMMON /DATA1/ IMAX.JKAX 
COMMON /DATA2/ XI , ERR, PI 
COMMON /DATA3/ ALP.APX 
COMMON /DATA4/ X(10) 

COMMON /DATAS/ XK1.XK2.DLT 
COMMON /DATAS/ A2.A3.A5.AG 



ft 


DO 100 1=1.10 
CALL DATA (ALPD) 

CALL NEWTON (Z.W) 

CALL VORTEX (W. GAM) 

C CALL PRESSURE (W, GAM) 

CALL LIST (GAM. W. Cl. 3) 

C CALL LIFT], (GAM . VI . CLW . CDW . CLF . CDF , CLT , CDT . CLL . 

C * TCL, TCD.RTOR. RTO) 

C WRITE (SC . 300) ALPD , CL3, CLW . CDW , CLF . CDF . CLT . CDT . CLL 
C * .TCL, TCD.RTOR. RTO 



I 


t ■' 
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ICO CONTINUE 

300 FORMAT ( IX , FS . 2 , IX , Dll . 4/5X , 6 ( IX . Dll . 4) /'5X , 5 ( IX . Dll . 4) ) 
200 FORMAT (2X. 'DLT' , 5X, ‘VORTEX POSITION* . 9X. 'GAM* , 9X. 

* ' CLL * , 9X , ' CLV ' , 9X . ' CLT ' /) 

STOP 

END 
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SUBROUTINE DATA (ALPD) 

INPUT OF DATA USED IN THE WHOLE PROGRAM 


IMAX 

JMAX 

XI 

ERR 

PI 

ALP 

APX 

X(I) 

XK 

DLT 


RNAK OF MATRIX i.e. THE NUMBER OF UNKNOWNS 
5 FOR NEWTON-PAPHSON METHODS (I-l.J) . (I+l.J) . (I . J*l) 

, (I.J+1) . (I.J) 

SQRT(-l) 

CONVERGENCE OF ROMBERG INTEGRATION 

3.14 

ANGLE OF ATTACK 

HALF OF APEX ANGLE OF MAIN DELTA KING 
INITIAL GUESS 

RATIO OF MAIN DELTA WING TO TOTAL SPAN LENGTH (=1) 
FALP DEFLECTION ANCLE 


A2 . A3 , A5 , AS : COEFF . OF SCHWARTZ - CHRI STOFFEL TRANSFORMATION 


IMPLICIT INTEGEP.*4(I-K) 

COMPLEX *16 XI. X 

REAL*3 ERR , P I , ALP . AFX , XKl , XK2 , DLT , A 2 , A3 , A5 , AS , ALPD , DLTD 

COMMON /DATA1/ I MAX , UMAX 

COMMON /DATA2/ XI, ERR, PI 

COMMON /DATA3/ ALP, APX 

COMMON /DATA4/ XflO) 

COMMON /DATA5/ XKl.m.DLT 
COMMON /DATAG/ A2, A3. AS. AS 
DATA ERR, PI / l.D-S, 3.141592654 / 

DATA XI / (0..1.) / 

DATA IMAX.JMAX / 1. 5 / 

DATA ALP, APX / .17453270 , .2685121925/ 

DATA XKl, XX2 / .6153946154. .8846152046 / 

DATA A2.A3.A5.A6 / -0.626892, -0.515958, 

- 0.978499. 1.3006 / 

WRITE (6.1) 

READ (5. *) DLTD 
READ (5, *) ALFD 
DLT=DLTD ‘PI/130. 

WRITE (6. 2) 

READ (5, " ) A2 . A3 , AS , A6 

WRITE { 24 . 100) ALP , APX . XKl . XK2 . DLT. A2 , A3 . AS . A6 
WRITE (6. 600) 

READ (5, *) (X(I) , 1=1, IMAX) 

FORMAT (IX, ’DLT*) 

FORMAT ( IX , ’ TRANS . COEFF . ' ) 

FORMAT (/5X, ' 1 . GIVEN DATA '/ 

' 6X, 'ANGLE OF ATTACK = '.Dll. 4/ • 

' 8X. 'APEX ANGLE = '.011.4/ 

' GX. 'FLPA RATIO = '.Dll. 4/ 

1 8X. 'TAP RATIO = '.Dll. 4/ 

' 8X, 'FLAP ANGLE = '.Dll. 4/ 

' 8X, 'TRAN. COEFF. : ' , 4 (2X.D11 . 4) /) 

45 
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600 FORMAT (IX, 1 INITIAL GUESS') 
RETURN 
END 


SUERCUTINE NEWTON (VZ.VVI) 


C 


2,00 


5 


6 

4 


CALCULATION OF VORTEX POSITION 

IMPLICIT INTEGER** (I -N) 

COMPLEX “16 XI , X.W, WN, Z, ZM, DEPv , DERM, VZ , VW 

REAL * 0 ERR .PI. ALP , AP X , XK1 , XK2 . DLT , A2 , A3 , A5 , A6 , SS , B , XZ , OA 

DIMENSION W(10.10) . 

- WN(10) ,2(10} ,DER{10) .ZM(iO.lO) .DERH(IO.IO) . 

' SS(IO.IO) ,B(10) , 

' IPVT(IO) ,0A(10,10) ,XZ(10) . 

■ VZ(10) .WJ(10) 

COMMON /DATA1/ IMAX.JMAX 
COMMON /DATA 2/ XI. ERR. PI 
COMMON /DATA3/ ALP.APX 
COMMO X /DATA*/ a (10) 

COMMON /DATA5/ XK1.XK2.DLT 
COMMON /DATAS/ A2.A3.A5.A6 
N=0 

CALL .TE'WPNTtX.W) 

DO 4 1=1. IMAX 
IX! 5 J-l.JMAX 
WN (J) =:W (I , J) 

CALL ZCALC (WN . 2 . DER) 

DO 6 J-l , JKAX 
ZM(I, J)=Z (J) 

DERM (I . J) =DEP. (J) 

CONTINUE 

CALL MAT EL (W . ZM . DERM . SS . B) 

CALL NORM (B , XLOG) 

IF(N.GT.SC) GO TO 11 
IF (XLOG . LT . -7 A CO TO 11 
CALL DECOMP (SS’ iPVT.OA) 

CALL SOLVE (SS.B.IPVT.XZ) 

CALL NSXFNT (X.XZ) 

WRITE (6 , 112) H. XLOG, (ZM(1. 5} , 1=1. IHAX) . (W(I.S) . 1=1. IMAX) 

N=N+1 
GO TO 100 

WRITE (6. 113) N, XLCC, (2M (I, 5) .1*1. IMAX) . (W(I.S) . 1*1. IMAX) 

DO 12 1*1. IMAX 
. VZ (I) -ZM (1 , 5) 

VW ( I ) =W ( 1 , 5) 

WRITE (34. 200) VZ(1).VW(1) 


41 * 


12 




113 


200 


FORMAT (IX. 12. 5X, ' MAXNRM = \D1S.G/ 

* 5X. ' Z= ' . 2D1S . 6/ 

* SX. 'W=\2D15.6) 

FORMAT (/5X. *2. VORTEX POSITION'. */ 

* ax.. 'VORTEX POSITION IN Z-PLANE 

* 8X. -’VORTEX POSITION IN W-PLANE 
RETURN 

END 


II 

ii! 

K j SUBROUTINE NEWPNT(XN.YN) 

L IMPLICIT INTEGER* 4 (I *N) 

I? ! COMPLEX* IS XN.YN. XI 

p j REAL*8 ERR.PI.XRE, XIM 

{’ j DIMENSION XN ( 10) . YN ( 10 . 10) 

| ! COMMON /DATA1/ XMAX.JKAX 

F | COMMON /DATA2/ XI . ERR . PI 

DO 1 1=1, 1 MAX 

I -; XRE=DREAL (XN (I) ) /100 . 

;/ XIM=DIMAG (XN ( I ) ) /100 . 

V YN(I.1)=XN(I)-XRE 

r i YN(I. 2)=XN(I) +XRE 

' j YN(I,3)=XN(I)-XI*XIM 

t YN ( 1 . 4) =XN ( I ) +XI *XIM 

• i 1 YN (1.5) =XN ( I ) 

> j RETURN 

l ! END 


1 

i 


'.2D11.4/ 
‘.2D11. V) 


Ki 

I-; 

J? 


r ■ 
r ■ 
( . 
K, 

t • 

i\ 


S 

I 

I s 

p 

& : 


SUBROUTINE MATED (X , Z , DER . S . BB) 

IMPLICIT INTEGER* 4 (I -N) 

COMPLEX* 16 X.Z. DER . XI . A. B , O, FUN2 

REAL * 3 S . BB . ERR . PI . ALP , APX . XK1 . XX2 , DLT . A2 . A3 . A3 , AS 

DIMENSION X(10.10) ,Z(10, 10) .DER (10. 10) ,3(10.10; .BB(IO) 
COMMON /DATA2/ XI. ERR. PI • 

COMMON /DATA 3/ ALP, APX 
COMMON /DATA 5/ XK1.XK2.DLT 
COMMON /DATA6 / A2.A3.A5.A6 
A= (FUN2 (X (1. 2) . Z (1 . 2) . DER ( 1 . 2) ) - 

* FUN2 (X(l, 1) .7.(1. 1) , DER (1,1)))/ 

* (X (1. 2) -X (1.1)) 

B= (FUN2 (X (1 . 4) . Z (1 . 4) . DER (1.4) ) •• 

* FUN2 fX (1, 3) ,Z(1, 3) .DER (1.3)))/ 

* ( (X (1,4) -X(l, 3) ) * (-XI) ) 

0=FUN2 (X(l.S) ,Z(1.5) .DER(l.S)) 
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S ( 1 , 1) =DREAL (A) 

S (1. -I) =DREAL (3) 

S ( 2 . 1) =DIMAG (A) 

S(2.2)=DIMAG(B) 

SB(1)=-DREAL(0) 

BB(2)=-DIMAG(0) 

RETURN 

END 


SUBROUTINE NORM (X . >ZLCCi ) 
IMPLICIT INTEGER *4(1 -H) 
ILLICIT REAL 4 8 (A-H.O-2) 
DIMENSION X(10) 

COMMON /DATA1/ IMAX.JWAX 

IM=IMAX*2 

XN0RM=O. 

DO 2 Jal.IM 

2 XNORM=XNORM + (DABS (X (J) ) } 4 * 2 

XL0G=DLCG10 (XNORM) 

RETURN 

END 


SUBROUTINE NEXPNT (X.Y) 

IMPLICIT INTEGER* 4 ( I -N) 

COMPLEX* 16 X.XI 
REAL *8 Y, ERR. PI 

DIMENSION a (10). Y (10) 

COMMON /DA.TA2/ XI. ERR. PI 
X(l) =X(1) +Y (1) *XX a Y (2) 

X ( 1) -IjREAL (X (1) ) -'•XI * (DABS (DIMAG (X (1) ) ) ) 

RETURN 

END 


SUBROUTINE ROMERO (IS. A, 3. RES) 

IMPLICIT INTEGER' 1 4 (I -N) 

COMPLEX * 16 A . D . RES . XI . Z . DEL , FUN1 . X . DI FF . FUNS 
REAL*8 ERR. PI 
DIMENSION 7 . (10. 10) 

COMMON /DATA 2/ XI. ERR. PI 
1=1 

DEL=B-A 

IF(IS.EQ.l) CO TO 111 

Z(1 . 1) = . 5* DEL* (FUNS (A) +FUN5 (B) ) 

GO TO 10 

111 2(1. 1) = .5*DEL* (FLINT (A) +FUN1 (B) ) 

10 J=IFIX (2 . * * (I- 1) ) 
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DEL=DEL/2 . 

1 = 1+1 

Z(I.1) = .5*Z(I-1.1) 

DO 1 K=l, J 

X=A+ (2 . • FLOAT (K) -1 . ) ‘DEL 
IF(IS.EQ.I) GO TO 3 
Z (1 . 1) =Z (1 . 1) +DEL*FUN5 (X) 

GO TO 1 

Z (1 . 1} ~Z (1 . 1) +DEL*FUN1 (X) 
CONTINUE 
DO 2 10*2,1 

CONTINUE 

DIFE"Z (I ,1) -Z(I.I-l) 
IF(CDA3S(DIFF) .LT.ERR) GO TO 20 
GO TO 10 
RES-Z(I.I) 

RETURN 

END 


FUNCTION FUN! (X) 

IMPLICIT INTEGER*-! (I-N) 

COMPLEX® 16 X. XI. FUN! 

REAL ' * S ERR . P I . XK1 . XK2 . DLT .A2.A3.A5.A6 

COMMON /DATA 2/ XI. ERR. PI 
COMMON /DATA5/ XKl.XK2.DLT 
COMMON /DATAS/ A2.A3.A5.AS 

FUN1=X* (X-AS) (DLT/PI) / (CDSQRT (X-A2) *CDSQRT(X-A6) 
• *(X-A3)** (DLT/PI)) 

RETURN 

END 


SUBROUTINE ZCALC (V. Z.DER) 

IMPLICIT INTEGER* 4 (I -N) 

COMPLEX a 16 V, Z.DER. XI , EPS, Sl'NC, B.DB. BP. BM. A2P.R3, RBP, 

* RBM . P2 . P 3P . P3M. XINT, XINT1 . XINT2 . XINT1P . XINT2P . 

* XINT1M . XINT2M . ZM . ZP . Cl . C2 . C3 
REAL*8 ERR.PI.XK1.XK2.DLT.A2.A3.A5.A6 
DIMENSION V{10) .2(10) .DER(IO) 

COMMON /DATA1/ IMAX.JMAX 

COMMON /DATA2/ XI.ERR.PI 

COMMON /DATA 5/ XX1.XK2.DLT 

COMMON /DATA6/ A2.A3.AS.A6 

EPS=0.01*XI 

C1=A2-A5 

C2=A2-A6 

C3=A2-A3 
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SING=2 . *A2*C1** (DLT/PI) "CDSQRT(EPS)/ 
* (CDSQRT (C?.) a C3* e (DLT/PX) ) 

DO 1 I=1.JHAX 
B=V (I) 

DB=B/100. 

BP=B*D3 

BM=B-DB 


A2P=A2*EPS 
RB=DREAL (B) 
RBP=DREAL (B?) 
RBM=DRi:AL (I’M) 
P2=A2*2 . *XI 
P3=RB*2.*XI 
P3P=RBP+2 . *XI 


* 

* 


P3M=RBM-*'2 . *XI 

CALL ROMBRG (1. A2P.P2. XINT) 
CALL ROMBRG ( 1 , P 2 .P3.XINT1) 
CALL ROMBRG (1.P3, B.XINT2) 
CALL ROMBRG ( 1 . P2 . P3P , XIKT1P ) 
CALL ROMBRG ( 1 . P 3P . B? . XINT2P ) 
CALL ROMBRG (1.P2.P3M,XINT1M) 
CALL ROMBRG ( 1 . F 3M , EM , XINT2M) 


a«5=XINrr*XINTlM+XINT2M+SING 
2 (I ) =XINT^XINT1 + XINX2+SING 
ZP=XINT*XINT1P+XINT2P*SING 

DER (I) « ( (CDLOG {(2(1) -ZM) / (2P-Z (I) )))/(?•• "03) ) “ 

(CDSQRT(8-AS) fl CDSQRT (3 -A2) a (B-A3) (DLi/Pl)/ 

(B«(B-AS)**(DLT/PI)>) 


CONTINUE 

RETURN 


END 


FUNCTION FUN2 (XI. Z.DER) 

IMPLICIT INTEGER* 4 (I -N) 

COMPLEX* 16 XI . Z , DER . XI , FUN 2 
REAL. ‘8 ERR. PI , ALP . APX , XK1 , XK2 , 

COMMON /DATA 2/ XI. ERR, PI 
COMMON /DATA 3/ ALP. APX 
COMMON /0ATA5/ XK1.XK2.DLT 
COMMON /BA7A6/ A2.A3.A5.A6 

, *DCONJG(Z) +XI * (XK1* (1 


DLT , A2 . A3 , AS , AS 


FUN2=APX # (: 


-(1 


-XK1) *CDEXP (-XI *DL7) ) ) / 
(XK1*ALP) 

1 /til /Xl-1 ./DCCNJG (XI) ) * (Xl-DCONJC (XI) ) ) ) 

; (CDS0RT{X1-A6) ‘CDSQRT (X1-A2) * (X1-A3) ** (DLT/PI)/ 
. (XI* (XI- AS) ** (DLT/PI) ) ) • 

* -DER/ (1 . /XI- 1 . /DCONJG (XI) ) 

RETURN 

END 
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SUBROUTINE VORTEX (W , GAM) 

IMPLICIT INTEGER* 4 (I -N) 

COMPLEX* 16 W.XI 

REAL *8 ERR , P I , GAM , ALP . APX 

DIMENSION W ( 10) , GAM (10) 

COMMON /DATA2/ XI, ERR. PI 
COMMON /DATA 3/ A LP.APX 

GAM(l) =2 . *PI •ALP/DREAL (XX * (1 ./W (1) -1 ./DCCNJG (W(i) ) ) ) 
WRITE (34. 100) GAM (1) 

ICO FORMAT (/5X, ' 3 . VORTEX STRENGTH'/ 

* 8X. 'VORTEX STRENGTH = '.Dll. 4/) 

RETURN 

END 


SUBROUTINE LIFT1 (GAM . VW . CLW, CDW, CLF , CDF , CLT, CDT, CLL 

* . TCL , TCD . RTOR , RTO) 

IMPLICIT INTEGER*4(I-N) 

COMPLEX* 16 VW , XI , X2 . X3 , X4, ES . A2P . A3M. A3? . ASM . 

* A5M. ASP. SINGS. SING6.SING7. SINGS. 

* C1.C2.C3.C4.C5.C6.C7.C8.C9.Z2.2. 

* XI . SIKG1 , SING2 . SING3. SING4. 21. Z2 . Z3, Z4. 

* CIO, Cll , C12, C13. C14.C1S, C16 , X5, X6, Z5. Z5 

REAL*3 GAM , A2 , A3 , AS , A6 . U , V , W , CP . ERR .PI, XK 2 , XK1 . DLT 

* ALP , APX , CL1 , 0.2 , CL3 . CL4 . CLW , CDW . CLF , CLT , CDT 

* CDF. CLL , RTOR . RTO , TCL , TCD . CLS . CL6 
DIMENSION VW (10) , GAM (10) , 2 (700) , CP (7C0) 

COMMON /DATA?./ XI. ERR, PI 

COMMON /DATA3/ ALP, APX 
COMMON /DAT AS/ XXI . XK2 , DLT 
COMMON- /DATA6/ A2.A3.AS.A6 
ES~ (A6-A2) /200 . 

A2P=A2+ES 

A3M=A3-ES 

A3P=A2*ES 

A6M-A6-ES 

A5M~A3~ES 

A5P=A5"ES 

C1=A2-A5 

C2=A2-A6 

C3-A2-A3 

C4-A3-A5 

C5=A3-A6 

C6-A3-A2 

C7=A6-A5 

C8-A6-A2 

C9=A6-A3 

C10=-A6 

CU=-A2 

C12=-A5 


£4 



C13=-A3 

C14=A5-A6 

C15»A5-A2 

C16*A5-A3 

SING1=2 . *A2*CDSQRT(ES) ‘Cl 46 (DLT/PI}/ 

* (CDSQRT (C2) "C3 6a (DLT/P I ) ) 

SING2=-PI ft A3‘C4“ (DLT/FI) 6 (-ES) 66 ( (PI-DLT)/PI)/ 

* ( (PI -DLT) ‘CDSQRT (C5) ‘CDSQRT (C6) ) 
SING3=PX*A3 6 C4** (DLT/PI) a E S 66 ( (PI -DLT) /PI) / 

* ( (PI -DLT) "CDSORT (C5) * CDSQRT (C6) ) 

SING4=-2 . ‘AS 6 C7“ (DLT/PI) ‘CDSQRT (-ES)/ 

6 (CDSQRT (C8) ®C9 * 4 (DLT/PI) ) 

SING5=-C12* C (DLT/PI) * (-ES) 4 *2/ 

4 (CDSQRT (CIO) ‘CDSQRT (Cl 1) ‘CIS 46 (DLT/PI) *2* 

SING6-C12** (DLT/PI) *ES**2/ 

4 (CDSQRT (CIO) "CDSQRT (Cll) *C13* 4 (DLT/PX) *2) 

SING7=-A5 4 PI*(-ES)* 6 ((PI«-DLT)/PI)/ 

* (CDSQRT (C14) ‘CDSQRT (CIS) *C1S‘ 4 (DLT/PI) 4 (PX+DLT)} 
SING8=A5 fi PI ‘ES* 6 ( (PI+DLT) /PI) / 

* (CDSQRT (C14) ‘CDSQRT (CIS) ‘CIS 4 " (DLT/PI) 6 (PI-DLT) ) 
DO 10 1=1,101 

X.l=A2P+ (A3M-A2P) *FL0AT(I-1)/100. 

CALL ROMBRG ( 1 , A2P , XI , Zl) 

Z(I)=Z1+SING1 

ZZ=Z(I) 

CALL VEL (W/ , CAM , XI . ZZ . U, V. VJ) 

10 CP(X)=1.a/U,P 6 ‘2-U 64 2-V‘ 4 2-W* 6 2 

X2*A3P* (- 1 357-A3P) ‘FLOAT (J. -102) /100. 

CALL ROMBRG ( 1 , A3P .X2.Z2) 

Z(I)=XI‘XK1*SING3*Z2 

ZZ=Z(I) 

CALL VEL (Wl , GAM . X2 . ZZ . U . V . Wt 
20 CP(I)=1.+ALP**2-U 4 *'2-V**2-H # *2 

DO 25 I 203. 303 

X3=- . 357+ (-ES+ . 357) "FLOAT (1-203) /ICO . 

CALL ROMBRG (1.A3P.X3.Z3) 

Z(I)=XI*XX1<-SINC3*Z3 

ZZ=Z(I) 

CALL VEL (VW.GAM.X3. ZZ.U. V. W) 

25 C?(I)=1.+A1P 4 ‘2-U 44 2-V* 4 2-W 44 2 

DO 30 1=304,404 

X4=ES+ ( . 4.75-ES) "FLOAT (1-304) /ICO. 

CALL ROMBRG ( 1 , ES , X4 , Z4) 

Z(I)=XI* (XK1+ (1 . -XK1) "CDEXP (XI ‘DLT) ) +SINC6+Z4 
ZZ=Z(I) 

GILL VEL (W.CAM.X4, ZZ.U, V.W) 

30 CP (I)=l . +AI.P* *2-’J“2-V“2-W 4c, 2 

DTi .5,0 T=AOC t:pc; 

XS= . 475* (ASM- . 475) *FLOAT(t-4O5)/lC0. 

CALL ROMERO (l.ES.XS.ZS) 


35 



40 


Z (I) =XI « (XK1+ (1 . -XK1) «CDEXP (XI *DL7) ) -:-SING6+25 
ZZ-Z(I) 

CALL VEL(VW.GAM.XS.ZZ.U.V.W) 

CP (I)=1.'>AL?* <> 2~U**2-V'* 6 2-H t< '2 
DO 50 1=505,606 

X6=A5P* (A6M-ASP) 'FLOAT (I -SOS) /ICO . 

CALL ROMRRG(l,A5P,X6,Z6) 

■ z(i)=xi*XKi+siNca+z6 
ZZ=Z(I) 

CALL VtL ( V VJ , GAM . X6 , ZZ , U . V , W) 

50 CP (I)=1.+ALP* < '2-U* 4 2-V**2-W®*2 

CALL SIM? (Z.CP.CLl, CL2,CL3,CL4,CL5,CL6) 

CLW= ( - CL6 •'•CL1 ) 

CDW-CLW® ALP 

CLF= (-CL5+CL2) *DCOS (DLT) * (1 .*ALP*APX e D7AN{DLT) ) 

CDF= (-CL5-S-CL2) *DCOS (DLT) * (ALP-A?X*DTAN(DLT) ) 

CLT= (-CL4+CL3) ‘DCOS (DLT) * (3. . +AL?*APX*DTAN (DLT) ) 

CDT= (-CL4+CL3) *DC0S (DLT) * (AL?-APX*DIAN (DLT) ) 

CLL=CLW+ (-CL5+CL2-CL4+CL3) *DCOS (DLT) 

TCL=CLW+CLF +CLT 
TCD=CDW+CDF+CDT 
RT0R=1./ALP 
RTO=TCL/TCD 
C WRITE (34, 100) 

C DO 46 1=1. -ICS, 5 

C46 WRITE (34, 2CQ) Z (I) , CP (I) 

va 1 TE { 3 4 , 300) CLW . CDW . CLF , CDF . CLT . COT , CLL . TCL , TDD . ESC , RTOR 
100 FORMAT (20X. ' Z * . 21X. ' Cp ‘ , 2X/) 

2 CO FORMAT ( 10X . 2D1 1 . 4 . 5X , D1 1 . 4) 

300 FORMAT {//IX. ' 6 . LIFT BY INTEGRATION OF PRESSURE'// 


* 

10X. 'LIFT OF MAIN WING 

= ',011.4/ 

* 

10-X, 'DRAG OF MAIN WING 

= '.Dll. 4/ 

* 

1CX. 'LIFT OF FLA? 

a I 
~ ’ « 

t 

10X. 'DRAG OF FLAP 

« » 

"" * 

# 

1CX. 'LIFT OF TAP 

= ’.Dll. 4/ 

* 

10X. 'DRAG OF TAP 

= '.Dll. 4/ 


10X, ' CHECK FOR LIFT 

= ’.Dll. 4/ 

* 

10X. 'TOTAL LIFT 

= '.Dll. 4/ 

4 

10X. 'TOTAL DRAG 

= '.Dll. 4/ 

* 

10X, 'LIFT DRAG RATIO 

-- ' .Dll. 4/ 

* 

10X. 'IDEAL RATIO 

= ' , Dll . 4) 


RETURN 


ENT) 


SUBROUTINE SIMP (Z, CP . CL! , CL2 . CL3. CL4. CL5, CL6) 
IMPLICIT INTECER*4(I-N) 

COMPLEX* 16 Z 

REAL'S CL1,CL2,CL3,CL4.C?,CL5,CL6 
DIMENSION Z (700) . CP (7C0) 





CL1=0. 

CL2=0. 

CL3=0. 

CL 4=0. 

CL5=0. 

CL6=0. 

do 10 1 = 1,100 

10 CL1=CL1 + CDABS (Z (1+1) -2 (1) ) 4 (CP (1*1) +C? (I) ) /2 . 
DO 20 1=102,201 

20 CL2=CL2+CDABS (Z (I+I) -Z (I) ) * (CP (1*1) +CP (X) ) /2 . 
DO 30 1=203,302 

30 CL 3=CL 3 + CDABS (Z (X+l) -Z (I) ) « (CP (1+1) +C? (I) ) /2 . 

DO 40 1=304,403 

40 CL4=CL4+CDA£S (Z (1+1) -Z (I) ) * (CP (1*1) +CP (I) )/2 . 
DO 50 1=405,504 

50 CL5=Cl.5+CDA3S (Z (1+1) -2(1) ) A (CP (I+l) +CP (I))/2. 
EX) 60 1=506,605 

60 CL6=CL6+CDA£S (2(1 + 1'.) -Z(I) ) * (CP (I+l) +CP (I) )/2 . 
RETURN 
END 


SUBROUTINE PRESSURE (VW, CAM) 

IMPLICIT INTEGER* 4 (I -N) 

COMPLEX* 16 VW, X. ES , A2P , ASM, A3? , ACM. Cl , 02 . C3, C4, C5. C6 . 

* C7.CS.C9.XI . 5 1 NCI , SING 2 ,'S IN03 , SXNC4 . Z1 . Z2 . Z 

REAL * 3 GAM , A2 . A3 , AS . A6 . U, V . W, C? , ERR . PI . XK1 . XK2 . DLT 

* .ALP.APX 
DIMENSION VW(10) ,GAM(10) 

COMMON /DATA2/ XI, ERR. PI 
COMMON /DATA3/ ALP.APX 
COMMON /DATA5/ XK1.XK2.DLT 
COMMON /DATA6/ A2,A3,A5,A6 
WRITE (34, 100) 

DO 200 1=1,41 

X=A2+ (A6-A2) 0 FLOAT (I -1) /40. 

IF (DREAL (X) . EQ.A2 . OR. DREAL (X) .EQ.A3.0R. 

* DREAL (X) . EO . O . .OR. DREAL (X) .EQ.A5.0R. 

* DREAL (X) .EQ.A6) CO TC 200 
ES= (A6-A2) /2C0. 

A2P=A2+ES 

A3M=A3-ES 

A3P=A3+ES 

A6M=A6-ES 

C1=A2-A5 

C2=A2-A6 


C3=A2 -A3 
C4-A3-A5 






, . \y-x- '*-. 


C5=A3-A6 

C6=A3-A2 

C7--A6-A5 

C8=A6-A2 

C9“AS-A3 

SING1=2>A2 4 CDSCRT(ES) ‘Cl* 4 (DLT/PI) / 
ft (CDSQRT (C2) ‘Ca*’* (DLT/PI) ) 

SING2 s =-PI Cl A3' 5 C4 4 ' > (DLT/PI) * (-ES) *“ ( (PI-DLT) /PI)/ 

® ( (PI-DLT) ftCDSQRT(CS) *CDSORT(C6) ) 

SING3=PI 4 A3*C4ft 4 (DLT/PI) ‘ES* 4 ( (PI-DLT) /PI) / 

4 ( (PI-DLT) ftCDSQRT (CS) *CDSQRT (OS) ) 

SINC4=-2. ftAS rt C7ft 4 (DLT/PI) *CDSQRT(-ES)/ 

4 (CDSQRT (C3) *C9ft 6 (DLT/PI) ) 

IF (DREAL (X) .LE .DREAL (A2P) ) GO TO 200 
IF (DREAL (X) . LT . DREAL (A3M) ) CO TO I 
IF (DREAL (X) . LE . DREAL (A3P) ) GO TO 2C0 
IF (DREAL (X) . I.T . DREAL (A6M) ) GO TO 2 
CO TO 200 

CALL ROMBEG ( 1 . A2P , X . Z1 ) 

Z*Z1+SIN51 
GO TO 300 

CALL RCMHRG ( 1 , A3P , X , Z2 ) 

Z--XI ftXKI+SING3-»-Z2 

CALL VEL (VW . GAM , X , Z , U, V, W) 

CF=1 . +ALP 4 * 2-U» * 2 - Vft o 2 -Hft 4 2 
WHITE (6. *) Z.C?,U,y,W 
WRITE ( 34, 4C0) Z . CP . U, V. W 
CONTINUE 

FORMAT {/5X, *4. PRESSURE DISTRIBUTION •/ 

* 14X. 'ZM0X, 'Cp'.lIX, 'U-M2X. 'V',12X, 1 W ' , IX) 

FORMAT (4X , 2D11 . 4 , 4 ( 2X, DU . 4) ) 

RETURN 

END 



t ■ 
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SUBROUTINE VEL(VVI.CAM,X.Z.U, V, W) 

IMPLICIT INTEGER * 4 ( I -N) 

COMPLEX* IS VW.A.FUN4.X.Z.XI 

REAL * 8 GAM .U.V.W, ERR, PI, ALP , A PX . XK1.XX2, DLT 

DIMENSION VW (10) .GAM (10) 

COMMON /DATA2/ XI. ERR, PI 
COMMON /DATA 3/ ALP.APX 
COMMON /DATA5/ XiCl.XK2.DLT 
A=FUN4 (VW ( 1 ) . GAM (1 ) . X) 

V^DREAL (A) 

WsDREAL (XI " A) 

IF (DREAL (X) . GT . 0 . ) GO TO 1 
U=1 . +DREAL (ALPftX+XI *GAH(1) * (CDLOG (X-VW(l) ) - 
ft CDLOG (X-DCCNJO (VW (1) ) ) -2. «PI »XI) / (2 . *PI) 

ft -DREAL ( Z ) 4 V-DIMAG f Z) *W) ‘APX/XK1 

CO TO 2 

1 U=4 . +DREAL (ALP 41 X vX 1 S GAM (1) * (CDLOG (X-VW ( 1) ) - 

se 



* 

* 


CDLCG (X-DCONJC (VW (1) ) ) ) / (2 . °PI) 
-DREAL (2) *V-DIMAG(Z) *W) *APX/XK1 

RETURN 

END 


FUNCTION FUN4 (V.G.X) 

IMPLICIT INTEGERM(I-N) 

COMPLEX* 16 V.X.FUN4.XI 

REAL*0 G .ERR , PI , ALP , APX . XK1 , XK2 , DLT , A2 , A3 , A5 , A6 

COMMON /DATA2/ XI. ERR. PI 

COMMON /DATA3/ ALP.APX 

COMMON /DATAS/ XXI , XK2 , DLT 

COITION /DATAS/ A2. A3, AS, AS 

FUN4= (ALP ■'■XI *G® (1 ./ (X-V) - 1 . / (X-DCONJG (V) ) ) / (2 . «PX) ) 

* *CDSQRT (X-A6) *CDSgRT (X-A2) • (X-A3) • • (DLT/PI) / 

* (X* (X-AS) (DLT/PI) ) 

RETURN 

END 


SUBROUTINE LIFT (GAM. W, CL) 

IMPLICIT INTEGER* 4 (I -N) 

COMPLEX » IS W . XI . A2P . A3M.A6M . ES . Cl , C2 . C3 . C4 , C5 . C5 . C7 . 

* C8, C9 , CIO, X1NTI , XINT2 . XINT3, SING1 , SING2 . 

* SING3.SING4.CLL 


REAL*Q 


GAM . CL 1 . CL2 . CL . ALP . APX , XK1 . XK2 . DLT . A2 . A3 . A5 . A6 , 


ERR. PI 


DIMENSION GAM(IO) ,W(10) 
COMMON /DATA 2/ XI, ERR. PI 
COMMON /DATA 3/ ALP.APX 
COMMON /DATAS/ XK1.XK2, DLT 
COMMON /DATA6/ A2.A3.AS.A6 
ES=(A6-A2)/100. 


C1=A2-A5 


C2=A2-A6 

C3=A2-A3 

C4=A3-AS 

CS=A3-A6 

C6=A3-A2 

C7-A6-A5 

C8=A6-A2 

C9=AG-A3 

C10=A5 

A2P=A2+ES 



. ■pnrgf.v 
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A3M=A3-ES 

ASM=A6-ES 

CLi=2 . ‘CAM (1) <>DIMAG (K (1) -DCONJG (H (1) ) ) «APX/ 

* (XKl 4 (XK1+ (1 . -XK1) ‘DCGS (DLT) ) ) 

CALL RCMBRG (2 , A2F , A3M.XINT1) ' ’ 

CALL EOMSRG ( 2 , A3P , CIO. XINT2) 

CALL KOMERG ( 2 , CIO , ASM, XINT3) 

SING1=2 . ‘A2‘ *2*01** (DLT/PI) ‘CD3QRT (ES) / 

* (CDSpRT (C2) *C3» * (DLT/PI) ) 

SING2--PI 4 A3* 4 2*04** (DLT/PI ) * (-ES) *• ((PI -DLT) /PI)/ 

*• ((PI -DLT) ®CDSQRT(C5) ‘CDSQRT(CG) ) 

SIUG3=-2 . ‘AS* (DLT/PI) *CDSQRT (-ES) / 

* (CDSQRT (CO) * C9 * 0 (DLT/? I)) 

S ING4=P X * A3 ‘ ® 2 ? C4 ® (DLT/PI) ( (PI -DLT) /PI) / 

* ( (PI -DLT) ‘CDSQRT (CS) ‘CDSQRT (C6) ) 
CLI^XIOT1+SING1*SING2*XINT3*SING3 

fc + (XIN72-5-SING4) ®DCOS (2 . ‘DLT) 

CL2=- 4 . 4 ALP ‘DIMAG (CLL) *APX/ (XXI* (XXI + (1 . -XK1) *DCOS (DLT) ) ) 
CL=CL1+CL2 

WRITE (6, *) CL1.CL2.CL 

WRITE ( 34 . 100) CL1 # CL2 , CL 

FORMAT ( /5X , * S . LIFT COEFFICIENTS . 1 / 


OX, • LIFT FROM VORTEX 


'.Dll. 4/ 


RETURN 

END 


SX, • LIFT FROM ATTACHED FLOW = ’,0X1.4/ 
QX. ’ TOTAL LIFT = ’.DU. 4 /) 


FUNCTION EUN5 (X) 

IMPLICIT INTEGER‘4(I-N) 

COMPLEX® 16 FUNS, X. XI 

REM, 1 3 ERR , PI , XXI , XK2 , DLT .hi. hi. MS , A6 

COMMON /DATA 2/ XI. ERR, PI 

COMMON /DhTAS/ XK1,XK2.DLT 

COMMON /DATA6/ A2.A3.A5.A6 

IF(DREAL(X) .EQ.O. .OR.DREAL(X) .EQ.A5) GO TO 1 

FUN5=X * * 2 4 (X-A5) “ (DLT/PI)/ 

* (CDSQRT (X-A6) ‘CDSQRT (X-A2) •> (X-A3) • * (DLT/PI) ) 

GO TO 2 
FuNS-0 . 

RETURN 

END 
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SUBROUTINE DECOM? (A. IPVT.OA) 

IMPLICIT INTEGERS f I -N) 

IMPLICIT REALMS (A-H.O-Z) 

DIMENSION A(IQ.IO) .WGRIC(IO) .OA(IO.IO) .2(10} ,IPVT (10) 
COMMON /DATA1/ IMAX. JMAX 
N-2 * IMAX 
DO 1 1=1. N 
DO 1 J=1,N 
GA(I.J)=A(I.J) 

CONTINUE 
IPVT (N) =1 

IF (N.EQ.l) GO TO £0 
NM1*N-1 

ANQRMsO. . 

DO 10 Jsl.N 
T=0. 

IX) 5 1=1. N 

T=T+DABS (A ( X , J) ) 

CONTINUE 

IF (T . GT . ANORM) ANORM=T 
CONTINUE 
DO 35 K=1.NM1 
KP1=K+1 

DO 15 I-KP1.N 

IF (DABS (A(I .K} ) .GT.DA3S (A(M.K) } } M=I 
CONTINUE 
IPVT(K)=M •••' 

IF (M - NE . K) IPVT (N) »- IPVT (N) 

T*A(M.K) 

A(M.K)=A(iC,K) 

A (K, K) =T 

IF (T.F.Q.O.) GO TO 33 
500 20 IssRPl.N. 

A(I ,K =-A(I.K)/T 
CONTINUE 
DO 30 J-KP1.N 
T=A(M,J) 

A(M, J) =A(K. J) 

A(K,J)=T 

IF {T . EQ . 0 . ) GO TO 30 
DO 25 I=KP1,N 

A(X,J)=A(I.J)+A(I.K) *T 

CONTINUE 

CONTINUE 
CONTINUE 
DO SO K=1.N 
T=0. 

IF (IC.EQ. 1) GO TO 4.5 
KM1=K-1 


Cl 
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40 

45 


50 


55 


DO 40 1=1. KM1 

T=T+A(I.K) 4 WORK (I) 
CONTINUE 
EK=1 . 

IF (T.LT.O. ) EK=- 1 . 
IF(A(K.IC) .EQ.O.) GO TO 90 
WORK (K) =- (EK+T) /A (K, K) 
CONTINUE 
DO 60 K3=1,NM1 
K=N-K3 
T=0. 

K?1=K*1 
DO 55 I=KP1.N 

T=T+A(I.K) 4 WORK (K) 
CONTINUE 
WORK(!C) =T 
M=IPVT (K) 


IF (M.EP.K} GO TO 60 
T=WORK(M) 

WORK (M) =WCRX (K) 

WORK(K)=T 
60 CONTINUE 
YN0RM=O . 

DO 65 1=1. M 

\7ICRM=YN0RM+DABS (WORK (I) ) 
65 CONTINUE 

CALL SOLVE (A, WORK, IPVT. Z) 
ZN0RN=0. 

DO 70 1=1, N 

ZN0RK=ZNGRM+DA3S (Z (I) ) 

70 CONTINUE 

CONG- ANCRM 4 ZNORM/YNCEM 
IF (COND .LT. 1 . ) C0ND=1. 
RETURN 
80 C0ND=1 . 

IF (A(l. 1) .NE.O.) RETURN 
90 C0ND=1 . E+32 
RETURN 
END 


SUBROUTINE SOLVE (A. B . IPVT. X) 

IMPLICIT INTEGER 4 4 (I -N) 

IMPLICIT P.EAL*8 (A-H. 0-Z) 

DIMENSION IPVT ( 10) , A (10. 1C) , 5 (10) , X (10) 
COMMON /DATA1/ IMAX.JMAX 
N=2* IMAX 
DO 1 1=1. N 
X(I)=3(I) 
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End of Document 



